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Abstract 

In this continuation paper the Schrodinger equation for the half-integral 
spin eigenfunctions is obtained and solved. We show that all the proper- 
ties already derived using the Heisemberg matrix calculation and Pauli's 
matrices are also obtained in the realm of these analytical functions. We 
also show that Einstein-Bose condensation for fermions is expected. We 
then conclude this series of two papers on the concept of classical spin. 

1 Introduction 

In paper VI (hereafter VI) of this series we developed a method by means of 
which it was possible to derive the classical coordinate-momentum representa- 
tion for the behavior of particles with half-integral spin. 

It was shown there that the matrix representation of this calculation is read- 
ily obtained when we pass from the active to the passive views and change the 
Poisson bracket to its similar matrix commutator. Then all matrix mechanics 
related to the spin property was derived by straightforward calculations. 

We will now make use of the fact that our classical calculations give us the 
coordinate-momentum representation of the spin property which is suitable to 
be quantized — in the Schrodinger or Dirac's sense — to derive the Schrodinger 
equation for the spin eigenfunction. This will be done in the second section of 
this paper. This will show that Heisemberg's matrix calculus and Schrodinger's 
analytical equation give identical results. 

In the third section the Schrodinger equation for the spin will be solved 
and we will show that all the properties derived by means of the algebraic 
calculations are also obtainable with analytical functions. 



1 



We will devote the fourth section to show quantitativelly that fcrmionic 
Einstein-Bose condensation is expected by the present theory. 
In the last section we will make our conclusions. 



2 The Spin Schrodinger Equation 



We are now in position to derive the Schrodinger equation that will enable us 
to obtain the half-integral spin eigenfunctions. 
The two functions 

S 2 and S3 (1) 

have to be written as operators in a quantization procedure. The spin eigen- 
function will be the function that makes these operators diagonal. 

The quantization procedure has to be undertaken with much care since the 
function S 2 has terms with products of position and momentum operators which 
do not commute. To see this one needs only to look at this function written in 
the coordinate-momentum representation 



S 2 = ^ 



(x 2 +y 2 ) 2 +P (p 2 x +P 2 y ) 2 + 2 (x 2 + y 2 ) [& + p 2 y ) 



This task is greatly simplified if we note that we might write 

„o 1 



(2) 



(3) 



where 



So — 




(4) 



This states the very difference from this problem to the one usually found in 
textbooks on the solution of orbital angular momentum Schrodinger equations. 
While in the later case one has only the operator L 2 with the operator L, its 
square root, unknown, in the present approach both are known by principle and 
their relation is defined by equation (^|). This means that we will have only to 
find the function that makes the operators So and S3 diagonal to automatically 
make S 2 also diagonal. 

The next step is to transform the coordinate-momentum representation of 
our problem into an operator representation. 

We begin with equation (El) and quantize the function So in the usual way, 



giving 



and also 




+pi)+r 



v 2 )\ 



(5) 



(6) 



2 



where 



(7) 



We might develop the product represented in expression (||) as a repeated 
application of a differential operator to get 



S 2 



I67' 



{[(Pi 



+ pl) 2 + 7 4 (x 2 + y 2 ) 2 + 2 7 2 (x 2 + y 2 ) (pi + f y ) 
-iikj 2 (xp x + yp y ) - Aj 2 h 2 } 



(8) 

where the crossed terms in momentum-coordinates were treated using Dirac's 
symmetrization procedure according to which 
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2 2 x /~-2-~2 1 ^2^2\ -~2-~2 o-fc^^ 

x P x -> 2 \ x Px + Vx% ) = x Px ~ 2thxp, 



(9) 



or simply developing the squared operator applied upon some function. The 
first four terms in expression (||) represent S 2 which means that we might write 
this expression as 



4 u 4 

as our final operator S 2 . 

Looking at expressions (uuj we see that, if equation 



Soip = HXip 



is satisfied, then the equation 

1 



S 2 ^ 



-81--)^ = * 



A 2 -l 



(10) 



(11) 



(12) 



is automatically satisfied and we have a relation between the eigenvalues of ([11]) 
and (H). 

Equation ([ll]) might be written, in a differential form and in rectangular 
coordinates, as 



^2 I P f & 2 _|_ ® 2 1 _j_ / L *( 2 -\- 1 2 \ ij) 
* a \dx 2 dy 2 ' ^ ; ' ' 



la 
1 



= h\ip. 



We might now introduce polar coordinates 

x = r cos 6 ; y = r sin ( 

and 

P 



(. 



1/4 



(13) 

(14) 
(15) 
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to find equation 



ia_ 

pdp 



dtp 



2A 



~p2~dT 2 



together with the equation for S3 given by 

— iU dip 



mip. 



2 d6 

Equation (|l7|) might be solved by putting 



t/j (p, 9) = e l2m0/R R{p) 



to get for equation ( |l6| ) 



ld_ 

pdp 



dR 
dp 



2A 



4m 



V> = o 



i? = 0. 



(16) 



(17) 



(18) 



(19) 



This equation will be solved in the next section. However, it is noteworthy 
that, for the wave function to be a continuous function of 9, when a rotation by 
an angle of 2kir radians is performed, we must have only half-integral values for 
m. 



3 Solution of the Equation 

The solution of equation ( |l9| ) follows a standard method which we will only 
briefly sketch here. We begin by looking at the asymptotic behavior of function 
R(p) for very large values of the variable p. For this case we put 



to find 



R(p) er k P 2 



-1/2 



as the only physically acceptable asymptotic behavior. 
For small values of the variable p we might try 



R(p) P -=5 p s 



to find 



s = ±2 \m\ ,±2 \m\ 
giving the four possible asymptotic behaviors 



(20) 
(21) 

(22) 
(23) 



i?(p)^P ±2|m| ,P ±2|mh \ 



(24) 
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where the first two choices select the even series by letting ao free (and making 
ai = 0) and the second two choices select the odd series by letting a\ free (and 
making ao = 0). We shall choose the plus sign on both cases because the wave 
function has to be finite at the origin p = 0. We then have 

R(p) p -^p*\™\, p 2\™\-\ (25) 

We next search for an expression for all values of p using the series expansion 

R(p)=p s e-^ 2 J2anP n - (26) 

n 

Substitution of expression (53) into equation ( |l9| ) gives the following equation 
5>" 2 { [(a + nf ~ 4 H 2 )] a nP n - 2 [A - (1 + s + n)] a n p n+2 ) = 0, (27) 

n 

which gives, for the coefficients, the recursion relation 

_ 2[(l + .s + n)-A] 

a n +2 - - — ; —2 7T-[2 a «- \ zii ) 

[n + s + 2) — 4 \m\ 

At this point it is easy to see that both choices of s in equation (^5|) will give 
exactly the same recursion relation and so, also the same series — direct substi- 
tution of s = 2\m\ with n even and substitution of s = 2 |m| — 1 with n odd 
shows this. We then will work only with the choice 

s = 2\m\. (29) 

If the series in expression (^6|) does not terminate, its asymptotic behavior 
when n —> oo is 

— (30) 

a n -2 n 

2 

which is the same as the term p n e p and is not acceptable as a physical asymp- 
totic behavior. Then the series shall terminate; this is accomplished by making 
the choice 

X N = (l + 2|m| +N), (31) 

for some value of n = N. Since N must be a positive number, we automatically 
find the relation 

H < ^y 1 - (32) 

The correct eigenfunction of the problem is given by 

JV 

R(p)=p 2 ^e-ip 2 Y / anP n (33) 

71 = 
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and is a solution of the system 



Comparing this last expression with 

S 2 ip = h 2 e{£+l)i;, (35) 

we find that 

t=^. (36) 

Equation ( |32"|) then means that 

|m| < £ (37) 

as expected. 

The multiplicity of our functions might be calculated with the use of the 
quantum number I. One might easily check that this multiplicity is given by 

2^+1 = Aat. (38) 

The final eigenfunction of our problem might be written in the (r, 9) repre- 
sentation as 

/ „ \ IH/ 2 - - ..,/. „ JL / „, \™/ 4 

- I __ I 2\m\ 2im9/h 



^ r,e)= \Ji?) r*We*™»i*e-*V*ir>» ) r 2^ (j^I J r n , (39) 

where the coefficients a n are given by expression ( p8| ) and the ratio a/ (3 is a 
structure constant used to identify the actual particle we are interested in — as 
seen in expression ([s]) . This structure constant is necessary since our calculations 
were general and reflect the behavior of any half-integer spin particle. It is then 
possible to calculate quantities such as the radius of the particle in terms of this 
constant by performing the integral 

f t ,\ m \ (a/ f3) = f 4 (r, 0)r^,M (r, 9)rdrd6. (40) 



As an example we might find the radius of a half spin particle (m = 1/2, 1 = 1/2) 
using its density function 

d(r) - \ip{r,9)\ 2 =iV 2 r 2 e -(*) 1/2 r \ (41) 

where iV is a normalization constant. This density is related with the internal 
structure distribution of the half-integral spin particles and shall not be inter- 
preted as probability distributions — we might get the mass distribution of the 
electron by multiplying the above expression by its mass, for example. 
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The mean radius of this particle is 

<ri/a - 1/a>= /°°r3e-»^r = — UJ ' ^ 

where we notice that the bigger the structure constant, the smallest the particle. 

The first possible values of A are given in table I. In this table we show the 
£ value related to the chosen A value and all possible values of |m|. The values 
of the cutoff number N and the multiplicity associated with I is also shown. 
It is remarkable that if we try to put odd values for variable A, or else, integer 
values for variable £, we cannot terminate the series in (|39|), since in this case the 
cutoff number N is odd, which is not allowed by our choice of the asymptotic 
behavior at the origin in expression ( p9| ) . Since the series does not terminate its 
asymptotic behavior for large values of variable p is no more given by expression 
( p0[ ) and we must reject these solutions. 

We also plot the density distributions for the cases [£ = 1/2, \m\ = 1/2), 
{£ = 7/2, |m| = 7/2) and (I = 9/2, |m| = 3/2) in figures 1, 2 and 3 respectively. 

Figure 1 shows that the £ — 1/2 spin particles structure is somewhat like a 
ring (zero density at the origin) with maximum density at a radius r e=1 / 2 .\ m \=i/2j 
depending on the structure constant. 

In figure 2 the same behavior of figure 1 is attained but we might see that 
the distance of the maximum density from the origin is now ^=7/2,^1=7/2(0/ P) 
which is greater than the one for the half spin case. We can visualize this as 
the increasing of a 'centrifugal' force giving the screening of the particle matter 
distribution — we are considering the same value for the structure constant. 

Figure 3 shows that the same ring structure will be present in all functions 
(they depend on r 2 '™' which always makes the densities tend to zero in the 
vicinities of the origin). For the cases where the difference in the quantum 
numbers £ and \m\ are different from zero we also find the appearance of nodes 
reflecting the multiple ring structure of these particles. The number of these 
concentric rings will be given by the expression (£ — \m\) + 1. 



4 Bose-Einstein Condensation 

In the previous paper (VI) we have shown that, if the parameter A has a lower 
bound, then the phenomenon of Bose-Einstein condensation will be expected 
for some temperature T con d- This is precisely what we have found. In Table I 
we note that the parameter A has the value A m i n = 2 as a lower bound. 
This implies that the fermion shall be supplied with at least the energy 

E min = 2hto (43) 

to continue to act as a fermion. 

Because of this lower bound for the energy of the fermions we expect that, 
at some value of the temperature, condensation takes place. 
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5 Conclusion 



In these two papers we aimed at showing that: (1) the concept of spin is not 
a particularity of the quantum mechanical formalism and might also be repre- 
sented in the realm of classical mechanics from where we can extract a model 
for it (or a picture). (2) even in this case 'space quantization' might be ob- 
tained apart from a constant (phenomenologically obtainable) which we relate 
to Planck's constant. (3) the concept of spin is related to the symmetries gener- 
ated by the Lie algebra associated with its Lie Group (SU(2)); since this group 
is the same generated by classical phase-space functions obeying the same Lie 
algebra with the product defined as the Poisson bracket, there is no impossibil- 
ity in deriving the concept of classical spin. (4) this classical representation of 
the spin might be 'quantized' using traditional methods to derive a Schrodinger 
equation which is an analytical representation for this quantity. The solution 
of this equation, when squared, will give us information on the particle internal 
structure, as for example its mass or charge distribution or its mean radius, in 
terms of some characteristic constant related to each particle and identifying 
it. (5) all expected quantum properties already obtained by Heisemberg matrix 
calculations, using Pauli's matrices, are also obtained with this method and 
this result reaffirms the formal identity between Hciscmbcrg's matrix calculus 
and Schrodinger's one. (6) spin is not a characteristic of relativistic calcula- 
tions although it is better represented in the realm of this theory where Lorentz 
invariance might be imposed. (7) fcrmionic Bose-Einstein condensation has a 
very intuitive explanation by means of the minimum energy the fermions have 
to possess in order to behave like fermions. 

We hope that these calculations will help in clarifying some misconceptions 
rather diffused in the literature about the classical versus quantum status of 
half-spin particles and also about the possibility of an analytical representation 
of the concept of half-spin. 

Aside the epistemological aspects, we also hope that the possession of the 
spin eigenfunctions will help the investigations in areas such as superconductiv- 
ity, Bose-Einstein condensation and many others. 
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Table 1: Values of A or I in terms of m and N. The multiplicity of each choice 
is also shown. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/quant-ph/9509013vl 
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